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We consider the nonspherically decaying radiation field that is generated by 

a polarization current with a superluminally rotating distribution pattern in 

— 1/2 

vacuum, a field that decays with the distance Rp from its source as R P , 
instead of Rp . It is shown (i) that the nonspherical decay of this emission 
remains in force at all distances from its source independently of the frequency 
of the radiation, (ii) that the part of the source that makes the main contri- 
bution toward the value of the nonspherically decaying field has a filamentary 
structure whose radial and azimuthal widths become narrower (as Rp 2 and 
Rp 3 , respectively), the farther the observer is from the source, (iii) that the 
loci on which the waves emanating from this filament interfere constructively 
delineate a radiation 'subbeam' that is nondiffracting in the polar direction, 
(iv) that the cross-sectional area of each nondiffracting subbeam increases 
as Rp, instead of Rp, so that the requirements of conservation of energy are 
met by the nonspherically decaying radiation automatically, and (v) that the 
overall radiation beam within which the field decays nonspherically consists, 
in general, of the incoherent superposition of such coherent nondiffracting 
subbeams. These findings are related to the recent construction and use of 
superluminal sources in the laboratory and numerical models of the emission 
from them. We also briefly discuss the relevance of these results to the giant 
pulses received from pulsars. © 2008 Optical Society of America 



1. Introduction 

LA. Preamble 

Maxwell's generalization of Ampere's law 1 establishes that electromagnetic radiation can 
be equally well generated by a time-dependent electric polarization current, with a density 
dP/dt, as by a current of accelerated free charges with the density j: 



here, E and H are the electric and magnetic fields, D is the displacement and c is the speed 
of light in vacuo. A remarkable aspect of the emission from such polarization currents is 
that the motion of the radiation source is not limited by c. Although the speed of charged 
particles cannot exceed c, nothing prevents the distribution pattern of a polarization current, 
created by the coordinated motion of subluminal particles, from moving faster than light. 2-4 
Indeed, radiation from such superluminal polarization currents has been observed in the 
laboratory. 5-8 
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Since electric polarization arises from separation of charges, a polarization current is by its 
nature volume-distributed. In fact, no superluminal source can be point-like; for, if a point 
source were to move faster than its own waves, it would generate caustics on which the field 
strength would diverge. 2 ' 9 

There is growing experimental and theoretical interest in radiation by polarization cur- 
rents whose distribution patterns move at a superluminal speed with acceleration. 8 One of 
the simplest implementations of such sources employs distribution patterns that have the 
time dependence of a traveling wave with circular superluminal motion; here, the acceleration 
is centripetal. We are investigating the use of polarization currents with such superluminally 
rotating distribution patterns in applications relating to communications and radar. 6,10 Fur- 
thermore, one of the proposed models of the radio emission from pulsars postulates the 
presence of sources of this type in the magnetospheres of rapidly rotating neutron stars. 11 ' 12 
The clarification of a diverse set of current questions, therefore, hinges on an understanding 
of the radiation from superluminal polarization currents undergoing circular motion. 13-15 

Our purposes in the present paper are (i) to examine the geometry of those regions within 
such extended sources that make the dominant contribution toward the radiation field ob- 
served at a given point and time, and (ii) to identify the salient features Uof the angular 
distribution of this radiation. A detailed knowledge of the extent and geometry of the con- 
tributing part of the source is required not only for the efficient design of practical superlumi- 
nal sources of this type (e.g. for the design of the dielectric in which the polarization current 
is generated), 6 but also for understanding the narrow widths of the giant pulses that are 
received from pulsars. 16 Likewise, a knowledge of the evolution of the angular distribution of 
the radiation with distance both facilitates the experimental detection of the tightly-beamed 
large-amplitude component of the emission from such sources and establishes a connection 
between two observed features (the nanostructure and the high brightness temperature) of 
the pulsar emission. 17-20 

In Ref. 13, the field of a superluminally rotating extended source was evaluated by su- 
perposing the fields of its constituent volume elements, i.e. by convolving its density with 
the familiar Lienard-Wiechert field of a rotating point source. This Lienard-Wiechert field is 
described by an expression essentially identical to that which is encountered in the analysis 
of synchrotron radiation, except that its value at any given observation time receives con- 
tributions from more than one retarded time. The multivalued nature of the retarded time 
is an important feature of all superluminal emission; we shall begin, therefore, by describing 
the relationship between observation (reception) time and retarded (emission) time for the 
particular case of a rotating source with the aid of Fig. [TJ 
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Fig. 1. (a) Cross-section of the Cerenkov-like envelope (bold curves) of the 
spherical Huygen's wave fronts (fine circles) emitted by a small element S 
within an extended, rotating superluminal source of angular velocity ui. S is 
on a circle of radius r = 2.5c/u, or, in our dimensionless units, f = ruo/c = 2.5; 
i.e. its instantaneous linear velocity is ruo = 2.5c. The cross-section is in the 
plane of S's rotation; dashed circles designate the light cylinder rp = c/uo 
(fp = 1) and the orbit of S. (b) Three-dimensional view of the light cylinder, 
the envelope of wave fronts emanating from S, and the cusp along which the 
two sheets <p± of this envelope meet tangentially. (c) The relationship between 
reception time tp and source (retarded) time t [ Eq. <^f] plotted for f = 2.5 
and three different observation points. The maxima and minima of curve (i) 
occur on the sheets <p± of the envelope, respectively. Curve (ii) corresponds to 
an observation point that is located on the cusp. Note that the waves emitted 
during an interval of retarded time centered at t c are received over a much 
shorter interval of observation time at tp c . Curve (iii) is for an observation point 
that is never crossed by the rotating sheets of the envelope. (After Ref. 13.) 
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l.B. Multivalued retarded times, the cusp and temporal focusing 

Parts (a) and (b) of Fig. [1] show the wave fronts that emanate from a small, circularly moving 
superluminal source S. As we have already pointed out, no super luminal source can be truly 
pointlike. Here we are considering a volume element of an extended source whose linear 
dimensions are much smaller than the other length scales of the problem. 

The emission of waves by any moving point source whose speed exceeds the wave speed 
is described by a Lienard-Wiechert field that has extended singularities. These singularities 
occur on the envelope of wave fronts where the Huygens wavelets emitted at differing re- 
tarded times interfere constructively and so form caustics. A well-understood example is the 
emission of acoustic waves by a point source that moves along a straight line with a constant 
supersonic speed. In this case, simple caustic forms along a cone issuing from the source, 
the so-called Mach cone, and most of the emitted energy is confined to the vicinity of this 
propagating 'shock' front. Another, similar example is the formation of the Cerenkov cone in 
the electromagnetic field of a uniformly moving point charge whose speed exceeds the speed 
of light inside a dielectric medium. 

When the supersonic or superluminal motion of such sources is in addition accelerated, 
the simple conical caustic that occurs in the Mach or Cerenkov radiation is replaced by a 
two-sheeted envelope with a cusp. 9,21 ' 22 The effect of acceleration is to give rise to a one- 
dimensional locus of observation points at which more than two simultaneously received 
wave fronts meet tangentially. The spherical wave fronts that are centered at the retarded 
positions of the source neighboring a point from which such coalescing wave fronts emanate 
cannot but be mutually tangential (in pairs) to two distinct surfaces, surfaces that constitute 
the separate sheets of a cusped envelope. 

More specifically, the Cerenkov-like envelope that is generated by a uniformly rotating su- 
perluminal source consists of a tube-like surface whose two sheets meet, and are tangent to 
one another, along a spiraling cusp curve; this envelope is depicted in Fig. [Hand mathemati- 
cally described in Eqs. flU])- f|TS]) below. At any given observation time, three wave fronts pass 
through an observation point inside the envelope, while only one wave front passes through 
a point outside this surface. The envelope and its cusp are the loci of observation points at 
which, respectively, two or three of the simultaneously received wave fronts are tangential 
to one another. To specify the retarded times t at which various wave fronts are emitted, 
let us adopt a cylindrical coordinate system based on the axis of rotation and denote the 
trajectory of the volume element S, shown in Fig. [1], by 

r = const., (f{t) — <p + out, z = const., (2) 

where (p denotes the initial value of ip, and uj is the angular velocity of S. Let a stationary 
observer be positioned at a point P, with cylindrical polar coordinates (rp, y?p, zp). The 
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retarded-time separation R(t) between the source volume element and the observer (i.e. 
their instantaneous separation at the time t of emission) will therefore be 

R(t) = [Op - zf + r\ + r 2 - 2rr P cos(y9 P - (p - wt)]* . (3) 

The relationship between the retarded time t and the observation time tp, i.e. 

c 

is plotted in Fig. [T](c) for the source speed ru = 2.5c and for three classes of stationary 
observation points: those, located sufficiently close to the plane of rotation, that are peri- 
odically crossed by the two sheets of the rigidly rotating envelope [curve (i)], or by just the 
cusp curve of the envelope [curve (ii)], and those at higher latitudes that are never crossed 
by the envelope [curve (hi)]. 

The ordinates of the neighboring extrema of curve (i) in Fig. [T](c) designate those observa- 
tion times, during each rotation period, at which the two sheets of the envelope go past the 
stationary observer [see Eqs. below]. Thus, the field inside the envelope receives con- 

tributions from three distinct values of the retarded time [curve (i)], while the field outside 
the envelope is influenced by only a single instant of emission time [curves (i) and (hi)]. The 
constructive interference of the emitted waves on the envelope (where two of the contributing 
retarded times coalesce) and on its cusp (where all three of the contributing retarded times 
coalesce [curve (ii)]) gives rise to the divergence of the Lienard-Wiechert field on these loci. 
There is a higher-order focusing of the waves, and so a higher-order mathematical singularity, 
on the cusp than on the envelope itself. While the singularity that occurs on the envelope 
is integrable, that which occurs on the cusp is not. In that it occurs in the temporal as 
well as the spatial domain, this focusing is distinct from that produced by a conventional 
horn, mirror or lens. The enhanced amplitude on the cusp is due to the contributions from 
emission over an extended period of source time reaching the observer over a significantly 
shorter period of observation time. 

The Lienard-Wiechert field derived in Ref. 13 was used as the Green's function for cal- 
culating the emission from a superluminal polarization current, comprising both poloidal 
and toroidal components, whose distribution pattern rotates (with an angular frequency uj) 
and oscillates (with a frequency Q) at the same time. 13 It was found that the convolution 
of the density of this current with the Green's function described above results in a field 
that decays nonspherically: a field whose strength diminishes with the distance Rp from the 

_i/9 1 

source as R P , rather than R P , within the bundle of cusps that emanate from the con- 
stituent volume elements of the source and extend into the far zone. This result, which has 
now been demonstrated experimentally, 6,7 was derived in Ref. 13 by setting the observation 
point within the bundle of generated cusps and evaluating the convolution integrals over 
various dimensions of the source. 13 The steps in this procedure are listed below. 
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1. The integration with respect to the azimuthal extent of the source was performed by 
means of Hadamard's method. 23 ' 24 It was shown that the Hadamard finite part of 
the divergent integral that describes the field of a superluminally rotating ring with a 
sinusoidal density distribution consists of two parts: one part is exclusively contributed 
by the two elements on the ring that approach the observer along the radiation direction 
with the speed of light at the retarded time (i.e. the elements for which dR/dt = — c), 
and the other part is contributed by the entire extent of the ring. 

2. The integration with respect to the radial dimension of the source was subsequently 
performed by the method of stationary phase. 25 

It was found that, when the radiation frequency is much higher than the rotation frequency 
uj, the main contribution toward the field of a superluminally rotating annular ring comes 
from the vicinity of the point on the ring that approaches the observer not only with the wave 
speed, but also with zero acceleration (i.e. the point at which dR/dt = —c and d 2 R/dt 2 = 
simultaneously). 

These contributing source elements are the ones for which the time-domain phase tp = 
t + R(t)/c is doubly stationary. Differentiating Eq. (j4j with respect to t, we can see that 

dR a d2R n r*\ 

— = -c and —0 (5) 

are equivalent to 

^ = and ^ = 0. (6) 

dt dt 2 v ; 

These conditions jointly define the point of inflection in curve (ii) of Fig. [U(c), corresponding 
to the cusp passing through the point of observation P. 

The collection of volume elements satisfying Eq. (J3j) within an extended source has a 
filamentary locus that is approximately parallel to the axis of rotation for an observation 
point located in the far zone (Fig. [2]). The nonspherically decaying field that is generated 
by a volume-distributed source arises almost exclusively from the elements in the vicinity of 
this narrow filament, a filament whose position within the source depends on the location of 
the observer. 

l.C. The zeroth-order evaluation of the angular position of the nonspherically decaying 
beam 

The angle of observation corresponding to the cusp, and the reason for the filamentary 
structure of the contributing parts of the extended source may be inferred from the above 
equations. Applying the first condition in Eq. (jSJ) to Eq. ([3]) and solving the resulting equation 
for the retarded time t, or equivalently the retarded position tp — out + 0, we obtain 

/l =f A^ 2 \ 

if = ip± = o9 P + 2tt — arccos — — , (7) 

V rr P / 
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Fig. 2. Schematic illustration of the light cylinder r = c/u, the filamentary 
part of the source that approaches the observeration point with the speed of 
light and zero acceleration at the retarded time, the orbit of this filamentary 
source, and the subbeam formed by the bundle of cusps that emanate from the 
constituent volume elements of this filament. The subbeam is diffractionless 
in the direction of Op. The figure represents a snapshot corresponding to a 
fixed value of the observation time tp. The polar width 50p of this subbeam 
decreases with the distance Rp in such a way that the thickness Rp50p of 
the subbeam in the polar direction remains constant: It equals the projection, 
Sz sin Op, of the z extent, 8z, of the contributing filamentary source onto a 
direction normal to the line of sight. The azimuthal width of the subbeam, on 
the other hand, is subject to diffraction as in any other radiation beam: 5<pp 
is independent of Rp. 



8 



where 

A = (f| -l)(r 2 -l)-(z-z P ) 2 . (8) 

In these expressions, (f, 5; fp, Zp) stand for (ruj/c, zu/c; rpu/c, zpu/c), i.e. for the coordinates 
(r, 2; rp, zp) of the source point and the observation point in units of the light-cylinder radius 
cjuj. (This radius, which automatically appears in the present calculations, turns out to be 
the main length scale of the problem.) 

The retarded times t± = (ip± — + 2n7r)/uj respectively represent the maxima and minima 
of curve (i) in Fig. [T^c) where n is an integer. Applying both conditions of Eq. (jHJ) to Eq. ([3]), 
we obtain Eq. ([7]) and A = 0. The retarded time t c = £±|a=o represents the inflection point 
of curve (ii) in Fig. [1(c). Curve (iii) in Fig. (T^c) corresponds to an observation point for 
which A < 0, and so if± are not real. 

The envelope of wave fronts comprises those observation points at which two retarded 
times coalesce, i.e. at which t = t±. Inserting these values of the retarded time in Eq. (jlj) 
and solving the resulting equation for fp as a function of (rp, zp) at a fixed observation time 
tp, we find that 

<fp = utp + f - 0±(r P , zp), (9) 

where 

/1 T AV 2 \ 

4>± = R± + 27r — arccos — — , (10) 

V rr P / 

with 

R ± = p - zpf + f 2 +fl- 2(1 T A 1 / 2 )] 1 ^. (ii) 

These equations describe a rigidly rotating surface in the space (rp, (fp, zp) of observation 
points that extends from the light cylinder fp = 1 to infinity (see Fig. CD). 

The two sheets <fi± of this envelope meet at a cusp. The cusp occurs along the curve 

A = 0, fp = utp + <p - <P±(rp,z P )\ A=0 , (12) 

shown in Fig. H^a). It can be easily seen that, for a far-field observation point with the 
spherical polar coordinates Rp = (rp + Zp) 1 ^ 2 , 6p = arccos(^p/i?p), fp, Eq. f|T2|) reduces to 

6p = arcsin(f~ 1 ) H , if P — if — §7r H , (13) 

to within the zeroth order in the small parameter -Rp 1 , where Rp = Rpu/c. [The higher 
order terms of this expansion are given in Eqs. (1551) and (|69p .] In other words, the cusp that 
is detected at an observation point (Rp, 9p,fp) in the far zone arises from the constructive 
interference of the waves that were emitted by the volume elements at f = esc 6p, ip = fp+^n, 
regardless of what their z coordinates may be. These volume elements therefore have a 
filamentary locus parallel to the axis of rotation whose length is of the order of the z extent 
of the source distribution along the line f = csc^p, if = fp + |7r (see Fig. EJ). 
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l.D. The filamentary locus of the contributing source elements 

The locus of source elements that approach the observer with the wave speed and zero 
acceleration at the retarded time has a filamentary shape not only within the zeroth order 
approximation in the small parameter -Rp 1 , but in general. To demonstrate this, we need to 
introduce the notion of bifurcation surface? 

When deriving the equation describing the envelope of wave fronts, we kept the coordinates 
(r,tp,z), which label a rotating source element, fixed and found the surface in the space 
(rp , (f-p , zp) of observation points on which dR/dt = — c at a given time tp. If we keep 
(rp, tp-p, zp) and tp fixed, then dR/dt = —c would describe a surface that resides in the space 
(r, (f, z) of source points: the so-called bifurcation surface of the observation point P. Like 
the envelope, the bifurcation surface consists of two sheets that meet tangentially along a 
cusp (a spiraling curve on which d 2 R/dt 2 = 0), but the bifurcation surface issues from the 
observation point P (rather than the source point S) and spirals about the rotation axis in the 
opposite direction to the envelope (see Fig. [3]). The similarity between the two surfaces stems 
from the following reciprocity properties of P and S: the equation describing the envelope, 
Eq. (Q, remains invariant under the interchanges r <-» rp , z «-> Zp , tp «-> — <pp, t <-» — tp. 

The locus of source elements that approach an observer P with the wave speed and zero 
acceleration at the retarded time is given by the intersection of the cusp curve of the bifurca- 
tion surface of P with the volume of the source. This filamentary locus has exactly the same 
shape as the cusp curve of the envelope [shown in Fig. H](a)], except that it resides in the 
space of source points, instead of the space of observation points, and points in the direction 
of the source velocity. The projection of this curve onto the (r, z) plane consists of a branch 
of a hyperbola with asymptotes that lie along the angles arcsin(f P " 1 ) and 7r — arcsin(fp 1 ) with 
respect to the z axis [see Fig. H(b)]. For an observation point that is located in the far zone, 
therefore, the projection of the cusp curve of the bifurcation surface onto the (r, z) plane is 
virtually parallel to the rotation axis. 

The reciprocity relations referred to above ensure that if a source element S is located on 
the cusp curve of the bifurcation surface of an observer P, then the envelope of the wave 
fronts emitted by S would have a cusp passing through P (or, conversely, if an observer P 
is located on the cusp curve of the envelope of wave fronts emitted by a source element 
S, then the cusp curve of the bifurcation surface of P would pass through S). In the case 
of a single point source, the retarded position tp of the source linearly changes with time 
(tp — tp + ut), and so the cusp that it generates is both spiral-shaped and rigidly rotates 
about the z axis. In the case of an extended source, on the other hand, the position tp of 
each contributing source element (an element that lies on the cusp curve of the bifurcation 
surface of a far-field observer P) is fixed (tp = tpp + 3ir/2, f = csc#p), and the elements that 
occupy that position are constantly changing. The cusps generated by the moving source 
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Fig. 3. The bifurcation surface of the observation point P for a source whose 
rotational motion is counterclockwise. The source points that lie inside this 
surface influence the field at P at three distinct values of the retarded time, 
while those that lie outside this surface influence the field at only a single 
value of the retarded time. The source elements on the filamentary locus at 
which the cusp curve of this surface intersects the source distribution approach 
P with the speed of light and zero acceleration at the retarded time and so 
generate a nonspherically decaying field at P. 
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elements that pass through this fixed position at various retarded times have a locus, at any 
given observation time, that is straight and stationary as shown in Fig. [21 In other words, 
the source elements constituting the filament at tp — ipp + 3ir/2, f = cscOp, each contribute a 
quasi-instantaneous 'pulse' of nonspherically decaying electromagnetic radiation that in the 
far field appears to have propagated out along a virtually straight-line locus defined by the 
angle Op = arcsin(f _1 ). 

I.E. Objectives and organization of the paper 

The objectives of the present paper are as follows (the location of the resolution of each 
objective is given in brackets): 

1. to show that the nonspherical decay of the radiation field that arises from a rotating 
superluminal source remains in force at all distances from this source independently of 
the frequency of the radiation (Section I3.DI) ; 

2. to specify the dimensions of the filamentary part of the source that makes the main 
contribution toward the value of the nonspherically decaying field [Eqs. (I59l)-(l60i)]; 

3. to show that the bundle of cusps emanating from this filament delineates a radiation 
'subbeam' that is nondiffr acting in one dimension; that is to say, the width of this 
beam in the polar direction remains the same at all distances from the source [Fig. [21 

Eq. my, 

4. to clarify how the requirements of conservation of energy are met by the nonspherically 
decaying radiation: the cross-sectional area of each nondiffracting subbeam increases 
as Rp, rather than Rp, with the distance Rp from the source (Section Hj); and 

5. to show that the overall radiation beam within which the field decays nonspherically 
consists, in general, of an incoherent superposition of the coherent nondiffracting sub- 
beams described above (Section HJ). 

We begin with the mathematical formulation of the problem in Section [21 In Section [31 
we show that objectives 1 and 2 can be achieved by replacing the method of stationary 
phase used in Ref. 13 with the method of steepest descents. 26 By converting the Fourier- 
type integral over the radial extent of the source to a Laplace-type integral and making 
use of contour integration, we present an asymptotic analysis for which the large parameter 
is the distance from the source (in units of the light-cylinder radius c/u) rather than the 
radiation frequency. Not only is there no restriction on the range of frequencies for which the 
emission from a rotating superluminal source decays nonspherically, but the more distant 
the observation point, the more accurate the asymptotic analysis that predicts this decay 
rate. 
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The more poweful asymptotic technique we employ here establishes, moreover, that the 
transverse dimensions of the filamentary part of the source responsible for the nonspherically 
decaying field are of the order of 5f oc Rp 2 in the radial direction and Sep oc Rp 3 in the 
azimuthal direction (see Section H]). The dimension of this filament in the direction parallel 
to the rotation axis is of the order of the length scale of the source distribution in that 
direction. 

The corresponding dimensions of the bundle of cusps that emanate from the contributing 
source elements can be easily inferred from the above dimensions of the filamentary region 
containing these elements. The cusps occupy a solid angle in the space of observation points 
whose azimuthal width 8(pp has a constant value (as does a conventional radiation beam) 
but whose polar width 50p decreases with the distance Rp as Rp 1 . This may be seen by 
considering a cohort of propagating polarization-current volume elements that are at the 
same azimuthal angle tp and radius r (possessing the same speed ru) but at differing heights 
z. Each will give rise to a cusp in the far zone that forms the angle 9p = arcsin(f _1 ) with the 
z axis, but starts from a different height at the light cylinder (see Fig. [2]). The spatial extent 
in the direction of increasing 9p of the composite set of cusps from this cohort of volume 
elements (the subbeam) will therefore be determined solely by the height 8z of the region 
confining the polarization current. Projected onto a direction normal to the line of sight, 
this will result in a width w = \5z\ sin^p occupied by the cusps that is independent of the 
distance Rp from the source. (Note that w is a fixed linear width, rather than an angular 
width.) 

Thus, the area Rp sin 0p50p<pp subtended by the bundle of cusps defining this subbeam 
increases as Rp, rather than Rp, with the distance Rp from the source. In order that the flux 
of energy remain the same across a cross section of the subbeam, therefore, it is essential 
that the Poynting vector associated with this radiation correspondingly decay as -Rp" 1 , rather 
than Rp 2 . This requirement is, of course, met automatically by the radiation that propagates 
along the nondiffracting subbeam. 

For a rotating superluminal source with the radial boundaries f < > 1 and f> > f<, the 
nonspherically decaying radiation is detectable in the far zone only within the conical shell 

arcsin(l/f > ) < Op < arcsin(l/f < ). (14) 

These limits on Op merely reflect the fact that a rigidly rotating extended source with finite 
radial spread entails a limited range of linear speeds ru; Eq. (Tl3l) shows that a limited range 
of speeds results in a limited spread in the angular positions of the generated subbeams. The 
overall beam described by Eq. (fl4|) consists, in general, of a superposition of nondiffracting 
subbeams with widely differing amplitudes and phases. The individual subbeams (which 
would be narrower and more distinguishable, the further away is the observer from the 
source) decay nonspherically, but the incoherence of their phase relationships ensures that 
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the integrated flux of energy associated with their superposition across this finite solid angle 
remains independent of Rp. 

Having made a preliminary description of the salient features of the analysis, we now 
embark on the detailed treatment of the problem in Sections [2] to HI We conclude in Section [5] 
with some remarks on the applicability of our analysis to numerical calculations of the 
emission from superluminal sources and to the observational data on the giant pulses received 
from pulsars. 

2. The nonspherically decaying component of the radiation field from a rotating 
superluminal source 

As in Ref. 13, we base our analysis on a polarization current density j = dP/dt for which 

Pr,tp,z(r, y?, z, t) = s rjtp>z (r, z) cos(m</?) cos(fit), -ix < <p < n, (15) 

with 

(p = if) — cut, (16) 

where P r ,ip, z are the components of the polarization P in a cylindrical coordinate system 
based on the axis of rotation, s(r, z) is an arbitrary vector that vanishes outside a finite 
region of the (r, z) space, and m is a positive integer. For a fixed value of t, the azimuthal 
dependence of the density ( TT5]) along each circle of radius r within the source is the same as 
that of a sinusoidal wave train with the wavelength 2irr/m whose m cycles fit around the 
circumference of the circle smoothly. As time elapses, this wave train both propagates around 
each circle with the velocity ru and oscillates in its amplitude with the frequency Q. This 
is a generic source: One can construct any distribution with a uniformly rotating pattern, 
P r ,ip,z( r , 0, z ), by the superposition over m of terms of the form s r}tp>z (r, z, m) cos{m<fi). 
The electromagnetic fields 

8A 

E = - v - 4 °-9(5ri- B = VpxA ' (17) 

that arise from such a source are given, in the absence of boundaries, by the following classical 
expression for the retarded four-potential: 

A"(x P ,tp) =c- 1 Jtfxdtf(^t)5(t v -t-R/c)/R yU = 0,---,3 (18) 

Here, (x P ,t P ) = (r P , ip P , zp, t P ) and (x, t) = (r,<p,z,i) are the space-time coordinates of 
the observation point and the source points, respectively, R stands for the magnitude of 
R = xp — x, and fi = 1, 2, 3 designate the spatial components, A and j, of A^ 1 and in a 
Cartesian coordinate system. 1 
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In Ref. 13, we first calculated the Lienard-Wiechert field that arises from a circularly mov- 
ing point source (representing a volume element of an extended source) with a superluminal 
speed ru > c, i. e. considered a generalization of the synchrotron radiation to the superlumi- 
nal regime. We then evaluated the integral representing the retarded field (rather than the 
retarded potential) of the extended source ()15p by superposing the fields generated by the 
constituent volume elements of this source, i.e. by using the generalization of the synchrotron 
field as the Green's function for the problem (see also Ref. 15). In the superluminal regime, 
this Green's function has extended singularities, singularities that arise from the constructive 
intereference of the emitted waves on the envelope of wave fronts and its cusp. 

Labeling each element of the extended source (1T51) by its Lagrangian coordinate (p and 
performing the integration with respect to t and if (or equivalently ip and if) in the multiple 
integral implied by Eqs. f|T5|) - f|T8l) . we showed in Ref. 13 that the resulting expression for the 
radiation field B (or E) consists of two parts: a part whose magnitude decays spherically, as 
Rp 1 , with the distance Rp from the source (as in any other conventional radiation field), and 

— 1/2 

another part B ns , with E ns = nxB ns whose magnitude decays as R P within the conical 
shell described by Eq. ( fT4|) . (Here, n = R/i? is a unit vector in the radiation direction.) 

The expression found in Ref. 13 [Eq. (47)] for the nonspherically decaying component of 
the field within this conical shell, in the far zone, is 

B- ~ -|iexp[i(fi/o;)(^P + 3tt/2)] £ m=/j± /iexp(-i/^ P ) 

x Y? j=1 q j f A > MrdzA- 1 / 2 u j exp(-if Jl (t)-), (19) 

(20) 



where 



H± = (n/u) ±m, (21) 
<fp = fp — utp, (22) 
qj = (l -iVL/uj in/u), (23) 



and 



s r cos ^pey + s^ej. 
Uj = { — s v cos^pey + s r e_|_ 
— s z sin 6*pei| 



(24) 



with j = 1, 2, 3. In the above expression, en = e^p (which is parallel to the plane of rotation) 
and = fix en comprise a pair of unit vectors normal to the radiation direction n. The 
domain of integration in Eq. (fl~9l) consists of the part of the source distribution s(r, z) that 
falls within A > (see Fig. HJ). 
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Fig. 4. (a) A segment of the cusp of the envelope of wave fronts emitted by 
a rotating point source with the speed ru = 3c. This curve is tangent to the 
light cylinder at the point (fp = 1, ipp = if — 3ir/2, zp = z) on the plane of the 
orbit and spirals outward into the far zone. Note that this figure represents 
a snapshot at a fixed value of the observation time tp. The cusp curve of the 
bifurcation surface of an observer P shown in Fig. [3] has precisely the same 
shape, except that it resides in the space of source points, instead of the space of 
observation points, and spirals in the counterclockwise direction: It is tangent 
to the light cylinder at the point (r = 1, ip — <fp + 37r/2, z = zp). (b) The 
projections of the cusp curve of the bifurcation surface and a localized source 
distribution onto the (f , z) plane. Only the part of the source that lies close to 
the cusp in A > contributes to the nonspherically decaying radiation. The 
source elements whose (f , z) coordinates fall in A < approach the observer 
with a speed dR/ dt < c at the retarded time and so make contributions toward 
the field that are no different from those made in the sub luminal regime. 
The asymptotes of the hyperbola A = make the angles arcsin(l/fp) and 
7r — arcsin(l/fp) with the z axis, so that for an observation point in the far 
zone (fp ^> 1) the projection of the cusp onto the (f, z) plane is (as depicted 
in Fig. [2]) effectively parallel to the rotation axis. 
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The function 0_(f , z) that appears in the phase of the integrand in Eq. (fT9l) is stationary 
as a function of f at 



f = f c (z) = {|(f| + 1) - [i(f| - l) 2 - (5 - 5 P ) 2 ] 1/2 } 1/2 - (25) 

When the observer is located in the far zone, this isolated stationary point coincides with 
the locus, 

f = f s = [l + (z-zp) 2 /(^-l)] 1/2 , (26) 

of source points that approach the observer with the speed of light and zero acceleration at 
the retarded time, i.e. with the projection of the cusp curve of the bifurcation surface onto 
the (r, z) plane (see Fig. H]). For Rp ^> 1, the separation f c — f s vanishes as Rp 2 [see Eq. fHUl) 
below] and both f c and f s assume the value csc^p. 
It follows from Eq. ([10]) that 



\r=r c 



= (f) c = R c + (p c - (p p , (27) 



_/df\? =f . c = 0, and 



<9f 2 



where 



and 



^a = -R^[{fl-l){fl -IT 1 -2], (28) 

f=f c 

fc = fp + 27r — arccos(fc/rp) (29) 



= r c (f 2 - rl) 1 ' 2 . (30) 

Note that for observation points of interest to us (the observation points located outside the 
plane of rotation, 6p ^ tt/2, in the far zone, Rp ^> 1), the parameter a has a value whose 
magnitude increases with increasing Rp: 

a ~ — Rp sin 4 # P sec 2 9 P (31) 

[see Eq. (1281) ]. In other words, the phase function 0_ is more peaked at its maximum, the 
farther the observation point is from the source. 

This property of the phase function <p_ distinguishes the asymptotic analysis that will 
be presented in the following section from those commonly encountered in radiation theory. 
What turns out to play the role of a large parameter in this asymptotic expansion is distance 
(Rp), not frequency (fi±ou). 
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3. Asymptotic analysis of the integral representing the field for large distance 

3. A. Transformation of the phase of the integrand into a canonical form 

The first step in the asymptotic analysis of the integral that appears in Eq. (fT9|) is to 
introduce a change of variable £ = £(f, z) that replaces the original phase <p_ of the integrand 
by as simple a polynomial as possible. This transformation should be one-to-one and should 
preserve the number and nature of the stationary points of the phase. 25 ' 26 Since 0_ has a 
single isolated stationary point at f = tq{z), it can be cast into a canonical form by means 
of the following transformation: 



r, z 



^ c (z) + la(z)e, (32) 



in which a is the coefficient given in Eqs. ( l28l) and ( 13T1) . 
The integral in Eq. ffTUl) can thus be written as 

/ fdfdzA- 1/2 u, exp(-i^_) = / dzd£ A(£, z) exp(i< 2 ), (33) 
ja>o h>is 

in which 

A(£,z) = fA" 1 / 2 u J |^exp(-i / i0 c ), (34) 

with 

^ = a£rR-(f*-l-A L ' 2 )- 1 , (35) 

and a = —fia/2. The stationary point f = tq and the boundary point f = rg respectively 
map onto £ = and 

e = es = _[ 2a -i(0 s _0 c )]V2 (36) 

where 

0s = 0_ | f=fg = 2tt - arccos[l/(f s f P )] + (f 2 r 2 - l) 1 / 2 . (37) 

The upper limit of integration in Eq. (!33l is determined by the image of the support of the 
source density (s in u.j) under the transformation (|32|) . 

The Jacobian dr/d£ of the above transformation is indeterminate at £ = 0. Its value at 
this critical point has to be found by repeated differentiation of Eq. (1321) with respect to £, 



df <9£ 



dr 

a£, (38) 



5 2 0_ ( df\ 2 86- d 2 f 



and the evaluation of the resulting relation ( 1391) at f = fc with the aid of Eq. (128j) . This 
procedure, which amounts to applying l'Hopital's rule, yields df/d^ = o — 1: a result that 
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could have been anticipated in light of the coincidence of transformation (|32|) with the Taylor 
expansion of 0_ about r — fc to within the leading order. Correspondingly, the amplitude 
A(^) that appears in Eq. (134")) has the value fc(fc — l) -1 u.j|^ = f» c exp(— i/x0c) a t the critical 
point C. 

When the observer is located in the far field (Rp ^> 1), the phase of the integrand on 
the right-hand side of Eq. (1551) is rapidly oscillating irrespective of how low the harmonic 
numbers fi± (i.e. the radiation frequencies /i±to>) may be. The leading contribution to the 
asymptotic value of integral (|33|) from the stationary point £ = can therefore be determined 
by the method of stationary phase. 25 However, in the limit Rp — > oo, £ s reduces to 

£ s ~ -3~ 1/2 cos 4 9p esc 5 OpRp 2 , (40) 

so that the stationary point £ = is separated from the boundary point £ = £s by an 
interval of the order of Rp 2 only. To determine the extent of the interval in f from which the 
dominant contribution toward the value of the radiation field arises, we therefore need to 
employ a more powerful technique for the asymptotic analysis of integral (1331) . a technique 
that is capable of handling the contributions from both fc and fg simultaneously. 

3.B. Contours of steepest descent 

The technique we shall employ for this purpose is the method of steepest descents. 26 We 
regard the variable of integration in 

I(z)= / €> d£A(£,5)exp(iaa (41) 

as complex, i.e. write £ = u + iv , and invoke Cauchy's integral theorem to deform the 
original path of integration into the contours of steepest descent that pass through each of 
the critical points £ = £s> £ = 0, and £ = £>. Here, we have introduced the real variable £>(£) 
to designate the image of f> under transformation fl32|) . i.e. the boundary of the support of 
the source term uj that appears in the amplitude A(£,z). We shall only treat the case in 
which fi (and hence a) is positive; I(z) for negative fi can then be obtained by taking the 
complex conjugate of the derived expression and replacing C with — C [see Eq. ( )34|) ]. 

The path of steepest descent through the stationary point C at which £ = is given, 
according to 

i£ 2 = -2uv + i(n 2 - v 2 ), (42) 
by u = v when a is positive. If we designate this path by C\ (see Fig. EJ), then 

p poo 

/ d£A(£, z) exp(i< 2 ) = (1 + i) / dvA\t =(1+i)v exp(-2av 2 ) 

J C\ J— oo 

~ (27r/ / u) 1/2 exp[-i(/i0 c - 7r/4)]uj| c csc0 P | sec ^(.Rp 172 , (43) 
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V 




Fig. 5. The integration contours in the complex plane £ = u + w. The critical 
point C lies at the origin, and u§ and u y are the images under transformation 
fl32|) of the radial boundaries f = fg(5) and f = f>(i) of the part of the source 
that lies within A > (see Fig. HJ). The contours C±, C2, and C3 are the paths 
of steepest descent through the stationary point C and the lower and upper 
boundaries of the integration domain, respectively. 



20 



for R-p ^> 1. Here, we have obtained the leading term in the asymptotic expansion of the 
above integral for large Rp by approximating A|^ = n + n„ by its value at C, where v = 0, and 
using Eq. (13T1) to replace a by its value in the far zone. Note that the next term in this 

* —1/2 

asymptotic expansion is by a factor of order R P smaller than this leading term. 

The path of steepest descent through the boundary point S, at which u = u$ = £s and 
v = [see Eqs. (1561) and pDj) ]. is given by u = — {v 2 + u s 2 ) 1//2 , i.e. by the contour designated 
as C*2 in Fig. [5j The real part of 

i£ 2 |c 2 = 2^ 2 + u s 2 ) 1/2 + m s 2 (44) 

is a monotonic function of v and so can be used as a curve parameter for contour C2 in place 
of v. If we let 2v(v 2 + u^ 2 ) 1 / 2 = —a, then it follows from 

£|c 2 = -(«s 2 + ^) 1/2 (45) 

that 

/ d£A(£,z)exp(i< 2 ) = !exp[i(cm s 2 -7r/2)] 
JC2 

x / o da(u s 2 + iay 1/2 A\^_ (us2+ia)1/2 exp(-aa). (46) 

The function A|c 2 that here enters the integrand can be determined only by inverting the 
original transformation ( |32f ). 

However, since the dominant contribution towards the asymptotic value of the above 
integral for Rp 3> 1 comes from the vicinity of the boundary point S, the required inversion 
of transformation (1321) needs to be carried out only to within the leading order in a. The 
Taylor expansions of 0±(r, z) about f = fs(i) are of the forms 

<t>± = 0s + rs^f 2 -l)(r 2 f 2 -l)- 1/2 (r-f s ) 

±§(2f s ) 3/2 (r 2 - l) 3 / 2 (f 2 f 2 - l)~ 3 / 2 (f - f s ) 3/2 + ■ • • . (47) 

According to Eqs. (1321) and (1361) . on the other hand, 

0_-0s = !a(£ 2 -£s 2 ). (48) 

In the vicinity of £ = £g, therefore, Eqs. (14TI) and (1481) jointly yield 

f ~ f s + § sin 5 # P sec 4 flpi^fe 2 - £ 2 ) (49) 

for Rp ^> 1. Note that £s 2 — £ 2 = — icr and that close to the cusp in the far zone 

A 1/2 ~ (2sin^p) 1/2 J R P (f -f s ) 1/2 , |f - f s | < 1, £ P » 1. (50) 
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Hence, inserting Eqs. (|49l) and (1501) in Eq. (1341) . taking the limit -Rp — > oo, and expressing £ 
in terms of a, we find that 

A\ C2 ~ exp[-i(/i0 c - 7r/4)]u i | s sin6'psec 2 6 l p(us 2 + icr) 1/2 o-" 1/2 (51) 

in the immediate vicinity of the point S at which a = 0. 

Strictly speaking, we should excise the singularity of A at o = by means of an arc-shaped 
contour. However, since this singularity is integrable and so has no associated residue, the 
contribution from such a contour vanishes in the limit that its arc length tends to zero. 
An alternative way of handling the removeable singularity at a = 0, followed below, is to 
introduce a change of integration variable. If we let a = r 2 , then the integral in Eq. fj46l) 
assumes the form 

/ d£A(£, z) exp(ia£ 2 ) ~ sin 6*p sec 2 #p exp[— i(/i0 c + 7r/4)]u,,| s / drexp(— ar 2 ) 
J Ci Jo 

~ i(27r/ / u) 1 / 2 csc^ P | sec# P | exp[-i(/^ c + tt^)]^-^^ 172 , (52) 

where use has been made of Eq. (13T1) and the definition a = — /ia/2. Note that this differs from 
the corresponding expression in Eq. (143]) for the integral over C\ by the factor | exp(— i7r/2) . 

The path of steepest descent through the boundary point £ = £>, at which u = «>, v = 0, 
is given by u = (v 2 + n > 2 ) 1 / 2 , i.e. by the contour designated as C3 in Fig. [5] The real part 
of the exponent 

i ? 2 |c 3 = -2^(^ + n > 2 ) 1 / 2 + in > 2 (53) 

is again a monotonic function of v and so can be used to parametrize contour C3 in place of 
v. If we let 2v(v 2 + n > 2 ) 1 / 2 = \i then it follows from 

e|c 3 = («> 2 + ix) 1/2 (54) 

that 

/ d£A(£, z) exp(iae 2 ) = \ exp[i(an> 2 - vr/2)] 

x J o dx(M> 2 + ix)" 1/2 ^l 5=(u> 2 +ix) i/ 2 exp(-«x)- (55) 

The asymptotic value of this integral for Rp 3> 1 receives its dominant contribution from 
X = 0. Because the function A\c 3 is regular, on the other hand, its value at x — can be 
found by simply evaluating the expression in Eq. ( l34l) at r = f>. The result, for i?p — > 00, is 

y4|c 3 , x =o — "r 2 sin 9p sec 2 6p(f 2 > sin 2 9p — l)~ Uj\f=f > exp(— ifx(pc)u > (56) 

[see Eqs. (jSJ)) and (I3ip]. This in conjunction with Watson's lemma therefore implies that 

f d£A(£, z) exp(i< 2 ) ^ 2 1/2 f 2 (f 2 sin 2 ^ P - l) _1 Uj| f> 

x exp[-i(/i0_| f> + 7r/ 4)])Lt — ^^p 1 , (57) 

to within the leading order in Rp 1 . 
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3. C. Asymptotic value of the radiation field 

The integral in Eq. (133!) equals the sum of the three contour integrals that appear in Eqs. (|43|) . 
(152]) and (l57j) ; the contributions of the contours that connect C\ and C*2, and C\ and C3, at 
infinity (see Fig. \5§ are exponentially small compared to those of Ci, C 2 and C 3 themselves. 
On the other hand, the leading term in the asymptotic value of the integral over C3 decreases 
(with increasing Rp) much faster than those of the integrals over C\ and Ci'- The integral 
over C 3 decays as Rp 1 , while the integrals over C\ and C 2 decay as Rp 1 ^ 2 [see Eqs. fH3|) . 
( 152]) . and ( 157]) ]. The leading term in the asymptotic expansion of the radiation field B ns for 
large Rp is therefore given, according to Eqs. ( 1T91) . ( 1331) . and ( 1431) . by 

B ns ~ -|(l + 2i)(27r) 1 / 2 J Rp 1/2 |sec0p|cscepexp[i(n/a;)(^p + 37r/2)]E j u= M± H 1/2 

xsgn(/i)exp(ifsgn/i)^ = igi/^ o d5u i | c exp[-i/i(0 c + &»)], (58) 

in which /x± can also be negative (see the first paragraph of Section 13. Bl) . 

This result agrees with that in Eq. (55) of Ref. 13. The two expressions differ by a factor 
of 2 — i because we have here included the additional contribution that arises from the source 
elements in the (vanishingly small) interval fs < f < tq. The integration with respect to 
f in Eq. (52) of Ref. 13 extends over fc < f < f>, while that in Eq. (14TI) extends over 
fs < f < The contribution from fs < f < fc is given, according to Cauchy's theorem, 
by the contribution from the lower half of C\ plus the contribution from C 2 . 

Even though the length of the interval f s < r < f c vanishes as i?p 2 as /2p tends to infinity 
[see Eq. (1401) ]. the contribution that arises from this interval towards the value of the field 

has the same order of magnitude as that which arises from fc < f < f>, and is by a factor 

~ 1/2 

of order R P greater than that which arises from the open interval fc < f < f>. Thus, 
the nonspherically decaying component of the radiation field that is observed at any given 
(xp, tp) arises from those elements of the source, located at the intersection of the cusp curve 
of the bifurcation surface with the volume of the source (Fig. HJ), that occupy the vanishingly 
small radial interval 

8f = fc — fs — I cos 4 #p esc 5 6pRp 2 (59) 

adjacent to the cusp at f = f g ~ csc^p [see Eqs. f l25l) and (126]) ]. 

The corresponding azimuthal extent of the source from which the contribution described 
by Eq. ( fl9|) arises is given by the separation + — 0_ of the two sheets of the bifurcation 
surface shown in Fig. [3] close to the cusp curve of this surface: The contribution of the source 
elements outside the bifurcation surface is by a factor of the order of Rp smaller than 
those of the elements close to the cusp inside this surface [see Eqs. (41) and (42) of Ref. 13]. 
Since <p + - 0_ ~ (2 5 / 2 /3)(csc # P )- 3 / 2 (f - f s ) 3/2 for |f - f s | < 1 and i? P > 1 [see Eq. (JEj)], 
and the contributing interval in f is of the order of Rp 2 [see Eq. ( 1591) ]. it follows that the 
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contributing interval in (p is 

8lf = (0+ - <f)-)r=r c ^ | cot 6 6 P Rp 3 . (60) 

The contribution from this vanishingly small azimuthal extent of the rotating source is made 
when the retarded position of this part of the source is (f = (fc [ see Eq. (I29p ]. i.e. when 
the contributing source elements approach the observer with the speed of light and zero 
acceleration along the radiation direction. Thus, the source that generates the nonspherically 
decaying field observed at a point (-Rp, Op, (fp) in the far zone (Rp ^> 1) consists entirely of the 
narrow filament parallel to the z axis that occupies a radial interval Sr oc Rp 2 encompassing 
f ~ csc^p and an azimuthal interval Stp oc Rp 3 encompassing if ~ Lfp + 3n/2 at the retarded 
time. 

3.D. Frequency independence of the nonspherical decay 

A further implication of the above analysis is that the generated field decays nonspherically 
irrespective of what the values of the frequencies fi±u may be. There is no approximation 
involved in introducing the transformation (132)) . and the asymptotic expansion is for large a. 
As derived here, therefore, the only condition for the validity of Eq. (158]) is that the absolute 
value of a ~ ^/i±Rp sin 4 6p sec 2 6p should be large, a condition that is automatically satisfied 
in the far zone for all nonzero frequencies. 

That the leading term in the asymptotic expansion of the integral in Eq. (133!) is pro- 
portional to Rp , instead of Rp , is a consequence of the particular features of the phase 
function 0_ described by Eqs. (f2~?j) - (l3"Ii) . These features originate in and reflect the particu- 
lar properties of the time-domain phase t + R(t) /c; they are totally independent of both the 
wavelength of the radiation and the size of the source. In contrast to all other nonspheri- 
cally decaying solutions of Maxwell's equations reported in the published literature (see, e.g., 
Refs. 27-29), whose slow spreading and decay only occur within the Fresnel distance from 
the source, the nonspherical decay that is discussed here remains in force at all distances. In 
fact, the greater the distance Rp from the source, the more the leading term dominates the 
asymptotic approximation in Eq. (158|) . 

The remaining z integration in the above expression for B ns amounts to a Fourier decom- 
position of the source densities s r ^ jZ \c with respect to z. Using Eqs. (l2TI) - (l30l) to replace <f>c 
in Eq. (158j) by its far-field value 

0c ~ £p - z cos Op + 3tt/2, (61) 

and using Eq. (1241) to write out Uj in terms of s r}V}Z , we find that the electric field E ns = fix B ns 
of the nonspherically decaying radiation is given by 

E ns ~ |(27r) 1 / 2 J Rp 1 / 2 |sec0p|csc0pexp[i(fi/a;)(^p + 37r/2)]E M = M± H 1/2 
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xsgn(/x) exp(i|sgnyu,) exp[-i/i(_R P + P + 37r/2)]{(is^ + Qs r /u;)e|| 

— [(is r — fis^/u;) cos^p + fis 2 sin#p/a>]ej_}, (62) 

in which s rj¥ , j2 stand for the following Fourier transforms of s r ip>z \c with respect to z: 

/oo 
dz s r ,<p,z(r, z) \, f=CSC Q p exp(i/i5 cos Op). (63) 
-oo 

This field is observable only at those polar angles Op within the interval arccos(l/f < ) < 
\0p — tt/2\ < arccos(l/f>) for which s r)iP)Z \r =csc e p are nonzero, i.e. at those observation points 
(outside the plane of rotation) the cusp curve of whose bifurcation surface (Fig. [3]) intersects 
the source distribution (Fig. Hj). 

3. E. Relevance to computational models of the emission from superluminal sources 

The asymptotic analysis outlined in this section provides a basis also for the computational 
treatment of the nonspherically decaying field B ns . The original formulation of B ns appearing 
in Eq. (119j) . in which the integral has a rapidly oscillating kernel, is not suitable for computing 
a field whose value in the radiation zone receives its main contribution from such small 
fractions of the f and (p integration domains as 5f oc Rp and 6<p oc -Rp 3 . The above 
conversion of the Fourier-type integrals to Laplace-type integrals renders the selecting out 
and handling of the contributions from integrands with such narrow supports numerically 
more feasible. 

4. The collection of nondiffracting subbeams delineating the overall distribution 
of the nonspherically decaying radiation 

We have seen that the wave fronts that emanate from a given volume element of a rotating 
superluminal source possess an envelope consisting of two sheets that meet along a cusp 
(Fig. [1]). There is a higher-order focusing involved in the generation of the cusp than in that 
of the envelope itself, so that the intensity of the radiation from an extended source attains 
its maximum on the bundle of cusps that are emitted by various source elements. If a source 
element approaches an observeration point P with the speed of light and zero acceleration 
along the radiation direction, then the cusp it generates passes through P. The reason is that 
both the locus of source elements that approach the observer with the speed of light and 
zero acceleration [i.e. the cusp curve of the bifurcation surface (Fig. [3])] and the cusp that is 
generated by a given source element are described by the same equation: The cusp curve of 
the bifurcation surface resides in the space of source points and so is given by Eq. (Tl2l) for 
a fixed (f P , ifp, zp), while the cusp curve of the envelope resides in the space of observation 
points and is given by Eq. (TT2l for a fixed (r, <p, z). 
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4- A. Angular extent of the nonspherically decaying emission 

The collection of cusp curves that are generated by the constituent volume elements of an 
extended source thus defines what might loosely be termed a 'radiation beam', although its 
characteristics are distinct from those of conventionally produced beams. The field decays 
nonspherically only along the bundle of cusp curves embodying this radiation beam. Since 
the cusp that is generated by a source element with the radial coordinate f lies on the cone 
Op = arcsin(l/f) in the far zone, the nonspherically decaying radiation that arises from a 
source distribution with the radial extent f< < f < r> is detectable only within the conical 
shell arcsin(l/f>) < Op < arcsin(l/f<). 

The field that is detected at a given point P within this conical shell arises almost exclu- 
sively from a filamentary part of the source parallel to the z axis whose radial and azimuthal 
extents are of the order of 5r oc Rp 2 and 5ip oc -Rp 3 , respectively [see Eqs. ( 1591) and (!60l) ]. 
The bundle of cusps emanating from this narrow filament occupies a much smaller solid 
angle than that described above. The parametric equation ip = zp(fp), tpp = <^p(fp) of the 
particular cusp curve that emanates from a given source element (r, ip, z) can be written, 
using Eq. ffT2l) . as 

zp = z±{rl -l)V2 (f 2_ 1) i/ 2; (64) 
(fp = ip — 2ir + arccos[l/(ffp)]. (65) 

If we rewrite Eqs. ( 1611) and (|65|) in terms of the spherical polar coordinates (Rp, 9p, ipp) of 
the observation point P and solve them for 9p and (fp as functions of (r, z) and Rp, we 
find that the cusp that is generated by a source point with the coordinates (r, (p, z) passes 
through the following two points on a sphere of radius -Rp: 

9p = arccos < — — 
[rRp 

and 

ipp = ip — 2n + arccos [l/(Rpf sin Op)}, (67) 

where the ± correspond to the two halves of this cusp curve above and below the plane of 
rotation (see Fig. H]). 

The Taylor expansion of the expressions on the right-hand sides of Eqs. ( )66|) and (IBTj) in 
powers of -Rp 1 yields 

Op = arcsin(l/f) - (z/^Rp 1 ± \{f 2 - If^Rp 2 + ■■■, (68) 

and 

(fP = <p - 3vr/2 - .Rp 1 + • ■ ■ . (69) 



% ± if 2 - l) 1/2 \Rl - 1 - t 



1/2 



(66) 
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These show that incremental changes 8r, 5z and 5(f in the position of the source element 
(r, (f>, z) result in the following changes in the (6p, ipp) coordinates of the point at which the 
cusp arising from that element intersects a sphere of radius -Rp in the far zone: 

S0 P = -f~\f 2 - l)' 1/2 5f - f-^zRp 1 + ■■■, (70) 

and Sipp = 5ip + ■ ■ ■. Because 5f is of the order of -Rp 2 , while 5z is of the order of unity for 
the filamentary source of the field that is detected at P, the dominant term in the expression 
on the right-hand side of Eq. (ITU]) is that proportional to the z extent of the filament. Given 
the observation point P, and hence a set of fixed values for the dimensions (5f, 5ip, 5z) of 
the filamentary source of the field that is detected at P, it therefore follows that the bundle 
of cusps generated by such a filament occupies a solid angle with the dimensions 

59p ~ — 5z sin OpRp 1 and 5(f P ~ 5tp (71) 

in the far zone. [Here, we have made use of the fact that 9p ~ arcsin(l/f) to within the 
zeroth order in Rp 1 to express f in Eq. (ITTj) in terms of 9p.] 

The bundle of cusp curves occupying the solid angle (ITTj) embodies a subbeam that does 
not diverge in the direction of 9p. The polar width S8p of this subbeam decreases with the 
distance Rp in such a way that the thickness RpSOp of the subbeam in the polar direction 
remains constant: It equals the projection of the z extent, Sz, of the contributing filamentary 
source onto a direction normal to the line of sight at all -Rp. The azimuthal width of the 
subbeam, on the other hand, diverges as does any other radiation beam: Sipp is independent 
of -Rp (see Fig. [2]). 

Thus, the bundle of cusps that emanates from the filamentary locus of the set of source 
elements responsible for the nonspherically decaying field at P intersects a large sphere of 
radius -Rp (enclosing the source) along a strip the thickness of whose narrow side is inde- 
pendent of .Rp. According to Eq. (ITT]) , the area Rp sin 9p59p5ipp subtended by this subbeam 
increases as -Rp, rather than Rp, with the radius of the sphere enclosing the source. Conser- 
vation of energy demands, therefore, that the Poynting vector associated with this radiation 
should correspondingly decrease as -Rp 1 instead of -Rp" 2 , in order that the flux of energy re- 
main the same across various cross sections of the subbeam. This requirement is, of course, 
automatically met by the (nonspherical) rate of decay of the intensity of the radiation that 
propagates along the subbeam. 

The nondiffr acting subbeams that are detected at two distinct observation points within 
the solid angle arcsin(l/f>) < Op < arcsin(l/f<) arise from two distinct filamentary parts of 
the source with essentially no common elements [see Eqs. (I59j) and (1601) ]. The subbeam that 
passes through an observation point P', though sharing the same general properties as that 
which passes through P, arises from those elements of the source, located at f = csc#p/, 
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tp' = (f-pr, that approach P', rather than P, with the speed of light and zero acceleration at 
the retarded time. Not only are the focused wave packets that embody the cusp constantly 
dispersed and reconstructed out of other (spherically spreading) waves, 9 but also the fila- 
ments that act as sources of these focused waves each occupy a vanishingly small (~ -Rp 5 ) 
disjoint volume of the overall source distribution, and so are essentially independent of one 
another. Unlike conventional radiation beams, which have fixed sources, the subbeam that 
passes through an observation point P arises from a source whose location and extent depend 
on P. 

It would be possible to identify the individual nondiffracting subbeams only in the case 
of a source whose length scale of spatial variations is comparable to Rp 2 (e.g. in the case 
of a turbulent plasma with a superluminally rotating macroscopic distribution). The overall 
beam within which the nonspherically decaying radiation is detectable would then consist 
of an incoherent superposition of coherent, nondiffracting subbeams with widely differing 
amplitudes and phases. The individual coherent subbeams decay nonspherically, but the 
incoherence of their phase relationships ensures that the integrated flux of energy associ- 
ated with their superposition across this finite solid angle remains independent of Rp. Note 
that the individual subbeams constituting the overall beam would be narrower and more 
distinguishable, the farther the observer is from the source. 

5. Concluding remarks 

The analysis we have presented here was motivated by questions encountered in the course of 
the design, construction, and testing of practical machines for investigating the emission from 
superluminal sources. 6 The original mathematical treatment of the nonspherically decaying 
radiation, 13 in which the integral over the volume of the source has a rapidly oscillating kernel, 
is not suitable for the computational modeling of the emission from such machines. We have 
seen that the nonspherically decaying radiation detected in the radiation zone receives its 
main contribution from such small fractions of the radial and azimuthal integration domains 
as Sr ~ Rp 2 and 5tp ~ Rp 3 . The above conversion of the Fourier-type integrals to Laplace- 
type integrals renders the selecting out and handling of the contributions from integrands 
with such narrow supports numerically more feasible. 

Not only the nonspherical decay of its intensity, but also the narrowness of both the beam 
into which it propagates and the region of the source from which it arises are features that are 
unique to the emission from a rotating superluminal source. These features are not shared by 
any other known emission mechanism. On the other hand, they are remarkably similar to the 
observed features of an emission that has long been known to radio astronomers: to the (as 
yet unexplained) extreme properties of the giant pulses that are received from pulsars (see, 
e.g. Refs. 16-20). The giant radio pulses from the Crab pulsar have a temporal structure of 
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the order of a nanosecond. 18 Under the assumption that they decay spherically like other 
conventional emissions, the observed values of these pulses' fluxes imply that their energy 
densities generally exceed the energy densities of both the magnetic field and the plasma 
in the magnetosphere of a pulsar by many orders of magnitude. 19 "The plasma structures 
responsible for these emissions must be smaller than one metre in size, making them by far 
the smallest objects ever detected and resolved outside the Solar System, and the brightest 
transient radio sources in the sky." 18 

The highly stable periodicity of the mean profiles of the observed pulses, 16 i.e. the rigidly 
rotating distribution of the radiation from pulsars, can only arise from a source whose dis- 
tribution pattern correspondingly rotates rigidly, a source whose average density depends on 
the azimuthal angle up in the combination (p — out only: Maxwell's equations demand that 
the charge and current densities that give rise to this radiation should have exactly the 
same symmetry (d/dt = —ood/dip) as that of the observed radiation fields E and B. On the 
other hand, the domain of applicability of such a symmetry casnnot be localized; a solution 
of Maxwell's equations that satisfies this symmetry applies either to the entire magneto- 
sphere or to a region whose boundary is an expanding wave front. Unless there is no plasma 
outside the light cylinder, therefore, the macroscopic distribution of electric current in the 
magnetosphere of a pulsar should have a superlumially rotating pattern in r > c/uo. The 
superluminal source described by Eq. (JT51) captures the essential features of the macroscopic 
charge-current distribution that is present in the magnetosphere of a pulsar and is thus an 
inevitable implication of the observational data on these objects. 

Once it is acknowledged that the source of the observed giant pulses should have a super- 
luminally rotating distribution pattern, the extreme values of their brightness temperature 
(~ 10 39 °K), temporal width (~ 1 ns), and source dimension (~ 1 m) are all explained by 
the results of the above analysis. The nonspherical decay of the resulting radiation would 
imply that the energy density and so the brightness temperature of the observed pulses are 
by a factor of the order of R/(r > — f<) 2 smaller than those that are normally estimated by 
using an inverse-square law, 13 a factor that ranges from 10 15 to 10 25 in the case of known 
pulsars. 16 

The nondiffr acting nature of this nonspherically decaying radiation [Eq. (1711) ]. together 
with its arising only from the filamentary part of the source that approaches the observer 
with the speed of light and zero acceleration [Eqs. ( |59l) and (J60l)]. likewise explain the values 
of its temporal width and source dimension. Furthermore, that the overall beam within 
which the nonspherically decaying radiation is detectable should in general consist of an 
incoherent superposition of coherent, nondiffracting subbeams (Section @J is consistent with 
the conclusion reached by Popov et al. 20 that "the radio emission of the Crab pulsar at the 
longitudes of the main pulse and interpulse consists entirely of giant pulses." 20 
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Two other features of the emission from a rotating superluminal source that were derived 
elsewhere 12 ' 30 are also consistent with the observational data on pulsars: 16 the occurrence of 
concurrent 'orthogonal' polarization modes with swinging position angles 12 and a broadband 
frequency spectrum. 30 
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